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O '. Abstract 



We present results for spectroscopic factors of the outermost shells in ^'^Ca 
and ^'^^Pb, which have been derived from the comparison between the avail- 



> 

o 
o 

' able quasielastic (e, e'p) data from NIKHEF-K and the corresponding calcu 



lated cross-sections obtained within a fully relativistic formalism. We include 
exactly the effect of Coulomb distortion on the electron wave functions and 



Q ' discuss its role in the extraction of the spectroscopic factors from experiment. 

^ ' Without any adjustable parameter, we find spectroscopic factors of about 

^ ■ 70%, consistent with theoretical predictions. We compare our results with 

' previous relativistic and nonrelativistic analyses of (e, e'p) data. In addition 

to Coulomb distortion effects we discuss different choices of the nucleon cur- 
rent operator and also analyze the effects due to the relativistic treatment of 

the outgoing-distorted and bound nucleon wave functions. 
25.30.-C, 24.10.Jv, 21.10.Jx 
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I. INTRODUCTION 



Spectroscopic factors and occupation probabilities are basic elements for our understand- 
ing of the nuclear structure, measuring the accuracy of the shell model description. The 
fundamental concept on which the shell model is based, the mean field approximation, is 
considered as the leading contribution in an expansion of multiparticle correlations within 
the general framework of the nuclear many-body theory. These correlations allow the nu- 
cleons to occupy partially different orbits above and below the Fermi level. Therefore, the 
deviation from full occupancy of the orbits below the Fermi level is a measure of the cor- 
relations neglected in the mean field approach, or in other words, a check of the validity of 
the mean field description. 

The spectroscopic factor Sa is defined as the probability to reach a final single particle 
(hole) state a when a nucleon is added to (removed from) the target nucleus. The occupation 
number A^^ is the number of nucleons in the quantum state a in the target nucleus, relative 
to the 2j + 1 limit. 

Theoretical spectroscopic factors and occupation numbers have been derived for some 
doubly magic nuclei, taking into account both short- and long-range correlations. In par- 
ticular, occupation numbers for the 3si/2 state in ^osp]^ have been obtained from different 
approaches. Mahaux and Sartor |I| derived an occupation number N^g = 0.83 and a spectro- 
scopic factor Sss = 0.69 by using a dispersion relation approach, which allows to extrapolate 
the empirical mean field from positive energies (optical potential) to negative energies (shell 
model potential). Slightly smaller spectroscopic factors have been obtained by Ma and 
Wambach 0] using a quasiparticle hamiltonian which includes correlations in a phenomeno- 
logical way. Pandharipande and collaborators obtained an average value 0.71 ±0.1 for 
the occupation probabilities of single-particle states just below the Fermi level, using a vari- 
ational calculation of nuclear matter to which a random phase approximation correction is 
added. 

From the experimental point of view, relative spectroscopic factors have been historically 
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determined by single-nucleon transfer reactions. In a search for absolute empirical values of 
spectroscopic factors, additional experimental information has been recently collected from 
two new sources. One is based on a sum rule analysis |^ of both transfer data and charge 
density differences of isotones extracted from (e, e') reactions. This method allows to go 
from relative spectroscopic factors to an absolute occupation number (see Ref. and Refs. 
therein). A typical value for the occupation number of the 3si/2 shell in ^''^Pb obtained from 
this method is Ns., = 0.78 ± 0.12 §. 

The other method is the quasielastic (e, e'p) reaction, which constitutes indeed a very well 
suited tool to extract experimental information on absolute spectroscopic factors. Modern 
facilities allow to study this reaction in detail and high precision measurements of cross 
sections are available [§]-§]. The small strength of the electromagnetic interaction allows 
one to study the transition in Born approximation (BA) by the exchange of a single virtual 
photon. Choosing properly the momentum and energy transferred by the virtual photon in 
order to satisfy the quasielastic condition ~ g^/2M), the process can be treated with 
confidence in Impulse Approximation (lA), i.e., assuming that the exchanged photon 
is absorbed by a single nucleon which is the one detected. Under these circumstances, it 
is possible to extract information on the energy and momentum distribution of the bound 
nucleon. 

Although the advantages of the (e, e'p) process over other possible reactions to study 
spectroscopic factors are widely recognized, the extraction of these factors from experiment 
is still not free of ambiguities. A reliable determination of spectroscopic factors requires 
an accurate knowledge of the mechanism of the reaction and, in this context, the exact 
treatment of the Coulomb distortion of the electrons is important, especially in heavy nuclei 
such as lead. 

Different methods have been proposed to handle the distortion. The early approaches 
were based on the eikonal approximation |[Ty] within a nonrelativistic scenario for the nuclear 
part. This method was applied to the analysis of (e, e'p) data taken at NIKHEF-K 
producing surprisingly low spectroscopic factors (0.49 ± 0.05 for the 3si/2 shell in ^°^Pb), 
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incompatible with results obtained from other reactions as well as with theoretical predic- 
tions. This fact, together with the large difference between the results from first and second 
order eikonal calculations, raised the question of whether this approximation was adequate 
to treat the Coulomb distortion or a more involved analysis was still needed. 

The first realistic calculations with a more exact treatment of the Coulomb distortion in 
{e,e'p) were made in 1990 by McDermott |[Tl|]. The incorporation of the distortion effects 
was achieved by means of a full partial wave analysis of the electron waves in the Coulomb 
potential of the target, as well as of the outgoing proton waves distorted by the optical 
potential. The complexity of the numerical calculations required a full relativistic framework, 
not only for the electron vertex where the energies involved in the experiments (hundreds of 
MeV) make it mandatory, but also for the nuclear vertex. Nevertheless, some approximations 
were still made (helicity conserved approximation, HCA) to simplify the treatment of the 
electron Coulomb distorted waves. The spectroscopic factor for the 3si/2 shell in ^ospj-^ 



amounts in this approximation to 0.65 manifestly above the value obtained within the 
eikonal approximation. 



Subsequently, a new relativistic calculation was reported by Jin et al. IT^, where the 



treatment of the Coulomb distortion of the electrons was exact and the above mentioned 
restricted approximation (HCA) was not made. The result of this calculation [|l2l for the 



spectroscopic factor in the same shell is 0.71, which is larger than the result of Ref. [|11 
The authors of Ref. [|1^ reported also some discrepancies with the results of Ref. [|11 
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the limit of plane waves for the electron (no Coulomb distortion) and stressed the need for 
further investigation. 

This paper is an attempt to clarify this situation by an independent analysis including 
the development of a new code |]T3|, which treats the Coulomb distortion of the electron in 
an exact way and uses a relativistic formalism for both the leptonic and the nuclear vertices. 
Our calculations are along similar lines to those of Ref. |12[ , though in some instances we use 



a different nucleon current operator and different bound nucleon wave functions. Differences 



and similarities with previous calculations in Refs. ||ll||r^ are discussed in detail in the next 



sections, which are organized as follows: In Sec. II we summarize briefly the formalism 
involved in our calculations. In Sec. Ill we discuss our results for the 3si/2 and 2^3/2 
shells in 208p|^ ^^^^ ^^le 2si/2 and 1^3/2 shells in ^°Ca. We discuss not only the effect 
of Coulomb distortion but also the effect of the relativistic optical potential and of nucleon 
current operator. Sec. IV summarizes the main conclusions. 



II. DESCRIPTION OF CALCULATIONS 

In this section we summarize the formalism used to describe the coincidence (e, e'p) 
reaction under the conditions defining the lA discussed in the Introduction. Fig. 1 represents 
graphically the (e, e'p) process. In this figure /cf (k^) is the four-momentum of the incoming 
(outgoing) electron and is the four-momentum of the exchanged photon. The ejected 
proton four-momentum is denoted by Pp. We use the notation and conventions of Ref. |T^ 



as well as h = c = 1. As represented in Fig. 1 the electromagnetic transition is treated 
in Born approximation, although we take into account the effect of the nuclear Coulomb 
potential to all orders by using electron distorted wave functions. 



The differential cross-section for this process is then written as |[13 

da _ S{ei + Ea - ef - Ep - Ea-i) 

defdVLfdEpdVLp ~ (2vr)5 



(2.1) 



if\ > 



where S indicates sum (average) over final (initial) polarizations and 

In this expression and stand for the electron and nuclear currents, respectively. The 
electron current is given by the well known point-like Dirac particle expression: 

jA'(r)=^^(r)7'^^Kr), (2.3) 
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where il^f,ipf stand for initial and final electron wave functions. In lA and within an inde- 
pendent particle model picture, the nuclear current can be written in terms of the nucleon 
current operator J'^ 



N 



J'Nir)=rArlJ'NrB{^, (2.4) 

with tp^, ipp the wave functions for the initial bound nucleon and final nucleon, respectively, 
and a nucleon current operator to be specified later. 

The initial and final electron wave functions, solutions of the Dirac equation with the 
Coulomb potential, have the form 



y.Y:^(h>)r.{f) , (2.5) 



where the functions 



(2.6) 



are eigenstates of total angular momentum with quantum numbers K/i (j = |k| — 1/2; Z = k 
if K > and / = — k — liffi;<0). The functions satisfy the usual radial equations 



15 1 and 4>'^{f) is given by 



r^f) = [Vi ® 5],^ = E < ^ ^ > yUr)xl . (2.7) 

For the outgoing electron wave function ip'j, the phase shifts 5^ have to be included with a 
minus sign. The functions fnidK and the phase shifts are obtained by numerical integration 
of the radial equations using the Milne procedure, as described in the work of Yennie et 
al. |jT6|, including up to third derivatives. The Coulomb potentials for ^osp]^ ^^^^^^ ^'^Ca are 
derived from the experimental charge distributions given in Ref. |]1^ . 

The bound state wave functions for the proton ip^ are spinors with well defined angular 



momentum quantum numbers kbI^b, and have a structure similar to that in Eq. 
They have been computed within the framework of the relativistic independent particle shell 



model. The mean field in the Dirac equation is determined through a Hartree procedure 
from a phenomenological relativistic lagrangian with scalar and vector S-V terms. We use 



the parameters of Ref . |]T8| , which are fitted to reproduce nuclear matter properties and the 
charge radius in ^°Ca, and the TIMORA code ||19[. 

The wave function of the detected proton ipp is a scattering solution of a Dirac-like 
equation, which includes S-V global optical potentials, obtained by fitting elastic proton 
scattering data [^. This wave function has basically the same structure of Eq. ( p. 51) except 



that, since the potential is in this case complex, the phase shifts and radial functions are also 
complex. In addition, since the wave function corresponds to an outgoing proton we have 
to use in Eq. ( |2.6| ) the complex conjugates /*, g* and in Eq. ( p.5|) the complex conjugate S* 
with a negative sign. 

For the current operator we consider the two choices ccl and cc2 introduced by de Forest 



21 in momentum space 



K W 

J^^{ccl) = {F, + RF,)r - ^{Pf + PiY , (2.8) 



Jj^(cc2)=Fi7^ + z^a^'^g., (2.9) 



where Fi and F2 are the nucleon form factors related in the usual way ||14[ to the electric 
and magnetic Sachs form factors of the dipole form. Pj in Eq. ( p.8|) is the four-momentum 
of the initial nucleon assuming on-shell kinematics [|21[ . 



As it is well known Eqs. (|2.8|) and (|2.9|) are equivalent when the initial and final 



nucleons are on-shell. For off-shell nucleons, as is our case, both expressions lead to different 
results and do not satisfy current conservation. In configuration space all three-momenta 
in Eqs. ( |2.8| ) and (|2.9| ) have to be considered as operators. For the cc2 choice, integration 
by parts allows one to replace the gradient operators acting on ?/^^, ipp by g*, the variable 
of integration in Eq. ( p.2| ). In the case of the ccl choice, it is not possible to get rid 
of the gradient operator acting on at least one of the nucleon wave functions. Given the 
fact that for off-shell nucleons none of the expressions (ccl, cc2) is fully satisfactory and 



both expressions fail to verify current conservation, when using the choice ccl we use the 
simphfying assumption of replacing Pp + Pj by the asymptotic Pp + Pm values, with P^ 
the missing momentum. This assumption simplifies enormously the calculations and is 
consistent with the prescription of de Forest in Ref. |2T| for the half-off- shell electron-proton 
cross-section a{ccl), which is commonly used. In the same spirit Fi and F2 are taken as the 
standard nucleon form factors at the asymptotic g-values, an approximation usually made 
when Coulomb distorted electron waves are employed ||TT| , [T2[| . Under this approximation the 
computation time is highly reduced. 



The operator cc2 was used in previous fully relativistic calculations [Tl|JT^. In this work 
we use both ccl and cc2 operators. Comparisons between half-off-shell electron-proton cross- 
sections obtained with both operators (ccl and cc2) have previously been made using plane 



waves for initial and final nucleons pl| , p2| . In Refs. and p2[ small differences (a few 



percent) were found at the kinematical conditions usually attained in the experiments. In 
the next section we present similar comparisons when realistic spinors are used to obtain 
the cross-section. 

Different notations have been introduced in the literature by different authors to distin- 
guish various degrees of approximations to the reaction mechanism. To avoid confusion we 
shall specify separately the approximation taken at each vertex, i.e., at the electronic and at 
the nuclear vertices. To distinguish the case when we use plane waves from the case when 
we use distorted waves in the electron current (Eq. (pl3|)), we use the notations PWBA and 
DWBA, respectively. On the other hand, to distinguish the cases when we use plane waves 
and distorted waves for the ejected proton, we use the notation PWIA and DWIA, respec- 
tively. Unless otherwise specified, the results presented are obtained with proton distorted 
waves (DWIA). In PWBA the electron current is given by 

jr,/.ee(r) = -^e^Ck^-'f^%kf, af)ru{h ^.) • (2.10) 

In the same way, when the optical potential is not included in the calculations (PWIA), the 
wave function of the ejected proton ipp becomes a plane wave in the corresponding equation 
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for the nuclear current. 

A few remarks concerning the numerical calculations in DWBA are in order. The trun- 
cation of the infinite sum over /t in Eq. ( |2.5| ) has been made to include at least 30 partial 
waves for both initial and final electron, as well as for the outgoing proton in DWIA. The 
radial integrals in Eq. (|2.2D have been carried out numerically up to typically 15 fm in the 



nuclear coordinate and 30 fm in the electron coordinate. This gives a good compromise 
in optimizing both numerical accuracy and computing time. The numerical accuracy has 
been checked by comparison to PWBA of the results obtained in DWBA for Z — > 0. The 



estimated numerical error amounts to less than 2% in the spectroscopic factors [13 



III. RESULTS 

We present our results in terms of reduced cross-sections p(-Pm) for selected values 
(i.e., for selected single-particle shells), defined by 

p{Pm) = I dE^ [a,p\PF\EFV , (3.1) 

JAE™ L -I dEpdefdilpdilf 

as functions of the missing momentum (P^ = Pa — Pa-i)- Experimentally, the integral over 
the missing energy {E^ = M^-i + M — M^) is taken over the interval Ai^^ that contains 
the peak of the transition under study. In our calculations of p{Pm), we take for crgp in Eq. 
( p.l|) the same expression used by the experimentalists, i.e., we use the expression acci given 
by Eq. (17) of Ref. [0]. 

We present results for the 3si/2 and 2C/3/2 shells in ^°^Pb and for the 2si/2 and 1^3/2 
shells in ''^Ca and compare with data in parallel kinematics (g || Pp)- We have chosen to 
study these shells for several reasons: i) they correspond to experimentally well separated 
peaks, ii) for these doubly magic nuclei the theoretical description is simpler and there are 
available optical potentials, in) this choice allows us to study different mass regions where 
the Coulomb distortion of the electron waves is expected to play a different role, and iv) 
for these shells we can compare our results with experimental data as well as with other 
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relativistic and/or nonrelativistic theoretical results. Unless otherwise specified the results 
presented correspond to parallel kinematics (g || Pp). All the calculations have been done 
for a fixed value of the kinetic energy of the outgoing proton (Tp = 100 MeV). We also take 
a fixed value of the incoming electron energy, = 412 MeV (375 MeV) for ^osp]^ (^°Ca). 



A. 3si/2 shell of ^'^^Pb 

We discuss results on the 3si/2 shell in detail because this shell is the most extensively 
studied in {e,e'p), both experimentally and theoretically. 

To start with, it is worth clarifying the situation with regards to the disagreement between 



the PWBA results in Refs. [|TI| and |T2|. As pointed out in Ref. [Il2|, both calculations 
being the same in PWBA should lead to the same results. We found that the reason for the 
disagreement is due to the fact that different wave functions were used for the bound proton. 
In our calculations we used the code by Horowitz and Serot (TIMORA) with the standard 
set of initial conditions and obtained a result (see dashed line in Fig. 2a) in agreement with 



that in Ref. [12]. With the same code but changing the initial conditions to those used by 



McDermott [ITll^il]) find a result (see dotted line in Fig. 2a) that agrees with the result 



in Ref. [|TT|. As indicated in the figure, the two above mentioned results have been obtained 



with the operator cc2, which is the operator always used in Refs. [pT] , p]2 . 

In what follows we use always for the bound proton the wave function obtained with the 
standard set of initial conditions. This wave function corresponds to a binding energy of 5.7 
MeV. 

In Fig. 2a we also show our result (solid line) in PWBA corresponding to the choice ccl 
for the current operator. As seen in this figure the PWBA results obtained with ccl (solid) 
and cc2 (dashed) current operators are very close. The main difference seen at the peaks is 
less than 8%, the ccl result being larger. 

In Fig. 2b we show our DWBA results obtained with the ccl operator (solid line) together 
with the DWBA results of Ref. |T2| (dashed line) that was obtained with the cc2 operator. 
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Also shown in this figure is the result obtained with the code of McDermott, used in Ref. 



Tl| , that makes use of the HCA (dotted line). In this figure, the three curves are obtained 
using the same relativistic optical potential |^ and also the same wave function for the 
bound proton. One can see that our results (solid line) and those from Ref. [|T2| (dashed 
line) are similar and show up the same qualitative changes with respect to the corresponding 
PWBA results shown in Fig. 2a. These changes are: i) a shift in Pm, ii) an increase at the 
maxima, and in) a filling of the minimum at Pm ~ 150 MeV. Note that the result obtained 
within HCA (dotted line in Fig. 2b) is not able to reproduce the two last mentioned effects 
(focussing effects) but it reproduces adequately the shift in Pm- 

In Fig. 3a we show in more detail the comparison between our PWBA and DWBA results 
shown by the dotted and solid lines, respectively. Also shown in Fig. 3a is the PWBA result 
with an effective g- value (dashed line), ge// = |^i,e// — ^/.e//!, where the electron initial and 
final effective momenta are given by 

7> r" n Za k , , 

fce// = A; + /. — ^ (3.2) 

with K = 1.1 A^^^ fm and fc = 1.35. This fc value has been adjusted to get the same 
Coulomb potential at the origin that corresponds to the charge distribution of the nuclear 
target obtained from the relativistic calculation, which agrees with the experimental one 



T7| . This value is somewhat smaller than the value fc = 3/2, corresponding to a spherical 
uniform charge distribution. In parallel kinematics, since Pm = Q — Pp is varied by varying 
g, a displacement in q produces a displacement in Pm- As seen in Fig. 3a, with this fc value 
the shift in Pm caused by electron Coulomb distortion is well accounted for. 

It should be pointed out that the effect of Coulomb distortion shows up differently in 
perpendicular kinematics (constant q and uj). To illustrate this point we show in Fig. 3b 
the comparison between our results in PWBA and DWBA in perpendicular kinematics. In 
this case there is no observable effect of displacement with Pm but there is a reduction of the 
cross-section in DWBA at Pm ~ that, as seen in the figure, is accounted for replacing q 
by qeff in PWBA. The reduction is due to the fact that the maximum of the bound nucleon 
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wave function (s-wave) at Pm = is not reached when q is replaced by Qeff- The actual effect 
of the focussing is of the same order in perpendicular kinematics than in parallel kinematics. 
On the other hand changing q into ge// produces quite a different effect in perpendicular 
kinematics because in this case Pm is varied keeping q constant and varying the angle between 
q and Pp. In what follows only the case of parallel kinematics is considered. 

The focussing effects, causing the filling of the minimum and increase of the maxima 
in Fig. 3a, are summarized in a quantitative way in Table I. In this table we quote the 
ratios between the reduced cross-sections calculated in DWBA and in PWBA at the two 
maxima. We compare the ratios obtained in this work with the ratios deduced from other 
relativistic and nonrelativistic calculations previously reported. The nonrelativistic results 
in Table I correspond to the code DWEEPY |10[. This code treats the Coulomb distortion 
in an approximate way making use of an expansion in powers of (Za/k) to first order 
(DWEEPYl) or to second order (DWEEPY2). DWEEPY2 is the most commonly used 
approximation in the nonrelativistic analyses of the data. One can see in Table I that our 
results and those of Ref. [O, which also treat Coulomb distortion in an exact way, are 



in very good agreement. The results obtained with the HCA code of McDermott ]TT|J2B 



do not reproduce adequately the focussing effect at the peaks. Among the nonrelativistic 
calculations, DWEEPYl overestimates largely the effect of focussing and DWEEPY2 gives 
still an overestimation. The spectroscopic factors are very sensitive to focussing effects 
because these factors are derived by scaling the theoretical reduced cross-section to the 
experimental one. As seen in Table I the various approaches to treat the electron Coulomb 
distortion fail to account for the focussing effect obtained with the exact DWBA calculation, 
which is non-negligible in ^°^Pb. Therefore, DWBA calculations are necessary to deduce 
reliable spectroscopic factors in heavy nuclei such as lead. 

In comparing relativistic and nonrelativistic calculations, it is also important to know 
how different approximations for the proton wave functions affect the reduced cross-section. 
To this end we compare in Fig. 4 the results of relativistic and nonrelativistic calculations 
in PWBA because in this case there is no effect from electron Coulomb distortion and the 
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differences come only from the various approximations at the nuclear vertex. 

In Fig. 4 the dotted and dashed lines show the PWIA results of the nonrelativistic and 
relativistic calculations, respectively. The solid and short-dashed lines show the results in 
DWIA of relativistic and nonrelativistic calculations. For the relativistic calculations we 



use the nuclear current operator Jf.c2 in Eq- ( P^-9D and the standard TIMORA solution |[n 
for the relativistic bound proton wave function. For the nonrelativistic bound proton wave 
function we use the upper component of the relativistic one, properly normalized, in the 
standard nonrelativistic form of the nucleon current, based on cc2 [^. The nonrelativistic 



optical potential has been taken from set II in Table 2.1 of Ref. 0, which was determined 



following the procedure of Ref. [^]. This set gives a good fit to elastic proton scattering 



data |26| from ^°^Pb at an energy of 98 MeV (which roughly corresponds to the proton 



energy involved in the {e,e'p) experiments 0). On the other hand, the relativistic optical 



potential used here was determined pO[ from a global fit to elastic proton scattering data 
from spherical nuclei with mass numbers 40 < A <208 in a wide range of projectile energies 
from 65 MeV to 1040 MeV. 

As seen in Fig. 4, both relativistic and nonrelativistic calculations give practically the 
same result in PWIA. The same is true when a Woods-Saxon potential is used for the bound 
nucleon [Q. This means that the use of nonrelativistic bound nucleon wave functions does 
not produce a significant change in the Pm range considered. 

The differences between relativistic and nonrelativistic calculations are much more no- 
ticeable in DWIA. As seen in Fig. 4, the shape is similar in both DWIA calculations but 
the relativistic potential leads to stronger absorption than the nonrelativistic one (about a 
20% difference at the peaks). This is at first sight surprising since both potentials reproduce 
well the elastic scattering of protons from this nucleus at proton energies of about 100 MeV. 
Therefore, the difference seen between the DWIA results in Fig. 4 must be attributed to 
details of the optical potential to which the elastic proton scattering at this energy is not 
sensitive. 

As it is well known [|l^ the relativistic S-V potential gives rise to non-local terms in the 
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nonrelativistic reduction. The effect of these terms in the (e, e'p) reaction was investigated 
by Boffi and collaborators 1^^. The authors of Ref. concluded that these terms only 



affect the inelastic processes, where the nuclear interior is important, producing about a 
15% increase of the absorption due to final state interactions in the (e, e'p) reaction. This 
conclusion is in agreement with the increased absorption observed in Fig. 4 in going from 
the nonrelativistic DWIA to the relativistic DWIA. It would be interesting to see to what 
extent the nonrelativistic reduction of the S-V optical potential may lead to a similar 



absorption as the relativistic DWIA. Work along these lines is in progress. 

B. Comparison with experiment 

The main results of this work are presented in Table II and in Fig. 5, where we compare 
our relativistic DWBA results to the experimental data, in parallel kinematics. Our theoret- 
ical results are given by the solid lines in Fig. 5 and have been scaled by the spectroscopic 
factors presented in Table II. 

Our spectroscopic factors for the 3si/2 and 2^3/2 shells in ^°^Pb, and for the 2si/2 and 
ld3/2 shells in ^°Ca, are given in the third and fourth columns of Table II, corresponding to 
calculations with ccl and cc2 operators, respectively. They have been obtained by scaling 
our theoretical results on p{Pm) to the experimental data from Refs. shown in Fig. 5. 
For each shell the overall scale factor (spectroscopic factor) has been obtained by means of 
an error weighted least squares procedure. The quoted errors in our spectroscopic factors 
include both statistical and systematical errors in experimental data Also given for 

comparison in Table II are the spectroscopic factors obtained from other relativistic and 



nonrelativistic calculations. The spectroscopic factors reported by Jin et al. have been 
obtained with a fully relativistic DWBA formalism similar to the one used in this work, 
using the cc2 current operator. As discussed in the previous section, in the case of the 



3si/2 shell in ^osp]^ authors of Ref. |]T2| use the same bound nucleon wave function used 



here. For ^°Ca the bound nucleon wave functions used in Ref. IT^ are different [E31 from the 
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standard TIMORA solutions used here. In column six of Table II we quote the spectroscopic 



factors obtained by McDermott ||TT| using also the cc2 operator. The last column of Table 
II contains the spectroscopic factors obtained with the nonrelativistic analyses reported in 
Refs. 1,0]. 

As seen in Table II, in general, the spectroscopic factors obtained with the cc2 current 
operator are somewhat larger than those obtained with cel. For ^''^Pb our cc2 results 
differ by less than 10% from the ccl results and are in agreement with the corresponding 
result by Jin et al. [|l^. As discussed in the previous section, the result by McDermott |TlI 



does not contain the correct focussing and corresponds to a different 3si/2 relativistic wave 
function leading to a smaller spectroscopic factor. Nevertheless, it is important to remark 
that all the results of the relativistic calculations for ^°^Pb are in agreement within a 10% 
and are compatible with theoretical predictions |l]-0], while the previous results from the 
nonrelativistic analyses are clearly too low. These too small spectroscopic factors in the 
nonrelativistic analyses of Refs. result from too much focussing of the electron waves 
and too little absorption of the outgoing proton wave, compared to the relativistic analyses, 
as was discussed in detail in the previous section (see Table I and Fig. 4). 

It should be stressed here that in our analysis the determination of the spectroscopic fac- 
tors has no free parameters. All of the parameters entering in the relativistic potentials were 
obtained from independent considerations [^,|l^]. Taking this into account it is remarkable 
the good quality of the fits to the experimental data seen in Fig. 5a for the two shells in 
^°^Pb. To have a measure for the quality of the fit, we consider the parameter Q defined 
as the x^-value divided by the degrees of freedom (number of data involved minus one) in 
the determination of the spectroscopic factor. In our cases Q < ?> corresponds to a good 
quality of the fit while Q > 5 corresponds to a poor fit. The Q-values obtained for ^osp]^ 
are Q < 3 for both shells (3si/2, 2^3/2) and for both types of current operators. Our results 
in Fig. 5a correspond to the calculations with the ccl operator, using the spectroscopic 
factors in column three of Table II. Similar results are obtained using the cc2 operator and 
the spectroscopic factors of column four in Table II. 
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In the case of ^"Ca the focussing effect is neghgible and the only noticeable effect of 
electron Coulomb distortion is a shift in P^- This is illustrated in Fig. 5b for the 2si/2 
orbital where we plot the DWBA result (solid line) and the result in PWBA with q^ff 
(dotted line), both using the ccl operator and the same spectroscopic factor {S = 0.44, see 
Table II). The effective momentum corresponds to fc = 1.35 in Eq. ( |3.2| ). As seen in the 
figure there is no significant difference between the results in DWBA and in PWBA with 
Qeff- Although not shown in the figure, the same is true with regards to DWBA and PWBA 
calculations for the 1^3/2 orbital in ^''Ca. 

As seen in Fig. 5b the agreement between our theoretical results for the 2si/2 and for 
the 1(^3/2 shells in ^°Ca with the ccl current operator is not as good as that obtained for the 
shells in ^'^'^Pb. If we use instead the cc2 operator, the agreement is not improved for the 
2si/2 shell, but it improves considerably for the 1(^3/2 shell. This is seen in Fig. 5b where we 
show our results for the 1*^3/2 shell with the ccl and the cc2 operators by solid and dashed 
lines respectively. For the case of the 2si/2 shell the Q value is larger than 7 in all cases, 
showing the poor quality of the fit. For the shell the Q value is almost 5 in the case of 
the ccl operator and is less than 2 in the case of the cc2 operator. 

It is important to realize that all the cases where our spectroscopic factors in columns 
three and four of Table II are low correspond to situations in which the shape of the the- 
oretical curve does not match well the trend in the experimental data and the Q values 
are large, i.e., the fits are poor. This is particularly the case for the 2si/2 shell in ^°Ca, 
where our spectroscopic factor is as low as that obtained in the nonrelativistic analysis of 
Ref. 1^. For the shell 1^3/2 in ^°Ca, where the fit with the cc2 operator is very good, we 
get a spectroscopic factor S = 0.76, similar to that found in ^''^Pb, and larger than the one 
obtained in the nonrelativistic analysis [|^. 

In the case of ^°Ca the comparison of our results with those of Ref. |J12|] is not as meaning- 
ful as in the case of ^ospb for the two following reasons: i) the bound nucleon wave functions 



used in Ref. []T2[ were different than the standard TIMORA solutions used here, and ii) 



the spectroscopic factors in Ref. |]12[ were derived by visual fitting. However it is interesting 
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to see that the result of Jin et al. [|T2| agrees with the one we obtain for the lc?3/2 shell with 
the cc2 operator, where we consider our deduced spectroscopic factor to be reliable. 



IV. CONCLUDING REMARKS 

We have analyzed the quasielastic (e, e'p) reaction within a fully relativistic formalism 
treating the Coulomb distortion of the electrons in an exact way. This analysis is expected to 
yield more reliable spectroscopic factors than previous analyses where the electron Coulomb 
distortion was not fully taken into account |TO|JTT[] . 



We present results for reduced cross-sections as functions of the missing momentum 
{Pm), corresponding to proton knock-out from the outermost shells in ^°Ca and in ^°^Pb. 
We study the effects purely due to electron Coulomb distortion as well as the effects due to 
the relativistic treatment of the initial (bound) and final (distorted) nucleon wave functions. 
We also study the effects of using different nucleon current operators. 

In parallel kinematics, Coulomb distortion produces two effects: i) displacement of the 
cross-section towards higher Pm values that can be simulated by the use of an effective 
momentum transfer ge// in PWBA and ii) focussing that shows up mainly in the maxima 
and minima of the reduced cross-section. 

The focussing effect plays an important role in the determination of the spectroscopic 
factors in ^°^Pb. In ^*^Ca the focussing effect is much smaller than in ^osp]^ j^^^g 
influence in the extraction of spectroscopic factors, but the displacement effect is still size- 
able. Exact DWBA calculations are worthwhile even in medium nuclei like ^°Ca in order to 
determine the precise magnitude of this displacement and to avoid introducing additional 
parameters. 

In perpendicular kinematics, we do not find displacement in Pm but there is an appar- 
ent increase of focussing (or antifocussing) at Pm ~ 0, caused by the effective momentum 
transfer. 

We find that the main difference between the results obtained with our relativistic treat- 
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ment of the nuclear vertex and with the usual nonrelativistic treatment [^,0 is due to the fact 
that the relativistic optical potential [^] produces more absorption in the (e, e'p) reaction. 
This gives rise to higher spectroscopic factors in much better agreement with theoretical 
predictions fTj-Q. In our view this is a strong point in favour of the relativistic S-V models 
that are able to account for both elastic proton-nucleus and (e, e'p) data simultaneously 
and confirms the fact that the relativistic phenomenology is in general superior to the usual 
nonrelativistic phenomenology in covering a large variety of experimental information on 
proton-nucleus scattering and particularly in what concerns spin rotation functions (see 



Ref. and references therein). 



We also find that the off-shell electron-proton cross-section is less sensitive to the choice 



of the nucleon current operators when using free Dirac spinors pT| , p2| than when using bound 
Dirac spinors. Thus, one should be cautious when estimating the uncertainty associated to 
the current operator from results based only on free spinors. Actually, further work is needed 
to clarify the situation with regards to the best choice of the nucleon current operator in 
the fully relativistic treatment of the (e, e'p) reaction. 

The spectroscopic factors have been obtained by scaling to experimental data the cal- 
culated reduced cross-sections without introducing any free parameter. For the 3si/2 and 
2d-i/2 shells in ^°'^Pb we find spectroscopic factors between 0.65 and 0.73 depending on the 
shell and on the nucleon current operator. High quality fits to the experimental reduced 
cross-sections are obtained in all the cases. These spectroscopic factors are in agreement 
with theoretical predictions and with the results obtained for the 3si/2 shell in ^'^^Pb by Jin 
et al. from a fully relativistic analysis analogous to the one carried out here. They are 



considerably larger than the spectroscopic factors obtained from the nonrelativistic analyses 
i- 

For '^'^Ca we only obtain a high quality fit to the data for the 1^3/2 shell when using the 
cc2 operator. In this case the value of the extracted spectroscopic factor is 0.76, similar to 
those obtained in ^°^Pb. For the 2si/2 shell the quality of the fit is very poor with both 
ccl and cc2 operators and, therefore, we cannot consider the deduced spectroscopic factors 



as being reliable. In our view, if the shape of the cross-section is not well reproduced, 
the derived spectroscopic factor is not meaningful. For this particular orbital the wave 
function obtained from the TIMORA code used here does not reproduce the experimentally 
observed shape of the reduced cross-section. In comparing with the results by Jin et al. 
for these two orbitals in ^°Ca, one should keep in mind that there are two main differences 
between their calculations and ours; i) the results reported in Ref. [|l^ were obtained using 
a phenomenological potential fitted to reproduce single-particle properties in ^°Ca ii) 
their spectroscopic factors are obtained by visual fit and no quality of the fit is reported to 
compare with ours. Nevertheless, in the case of the shell 1^3/2 with the cc2 operator, where 



the quality of our fit is good, our spectroscopic factor agrees with that of Ref. |12 



In this work we made no attempt to fit to experiment the r.m.s. radii and binding energies 
of the orbitals. Thus, the fact that we get high quality fits for all the orbitals studied, except 
the 2si/2 in "^^Ca, can be considered as a success of the relativistic analysis. It would be 
interesting to explore whether more elaborated relativistic S-V models, including non linear 
terms, may produce results in better agreement with experimental data on reduced cross- 
1/2 



sections for the 2si 10 shell in ^''Ca. 
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FIGURES 



FIG. 1. Schematic picture of the {e,e'p) process. 

FIG. 2. Reduced cross-sections for the shell 3si/2 of ^''^Pb in DWIA. (a) Comparison of PWBA 
results with ccl and cc2 operators. The dashed and dotted lines correspond to different choices of 
the 3si/2 wave function, (b) Comparison of DWBA results with ccl and cc2 operators. See text. 

FIG. 3. Reduced cross-sections for the shell 3si/2 of ^^^Pb in DWIA corresponding to parallel 
(a) and perpendicular (b) kinematics. The results obtained in DWBA are compared to the results 
in PWBA with and without effective momentum transfer. 

FIG. 4. Comparison of PWBA results from relativistic (r.) and nonrelativistic (n.r.) calcula- 
tions in DWIA and PWIA for the 3si/2 sheU in ^ospb. 

FIG. 5. Theoretical and experimental reduced cross-sections for the 3s 1/2 and 2^3/2 shells 
in ^'^^Pb (a) and for the 2si/2 and 1^3/2 shells in ^''Ca (b). Experimental data are from Refs. 
I^Q. Solid lines are the results of DWBA+DWIA calculations with ccl operator, scaled by the 
spectroscopic factors in the third column of Table II. The dotted line in (b) corresponds to the 
result with PWBA+DWIA with the ccl operator. The dashed line in (b) corresponds to the result 
with the cc2 operator. 
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TABLES 

TABLE L Ratio between DWBA and PWBA reduced cross-sections at the two maxima in the 
3si/2 shell of 208pb ( Tp = 100 MeV and ej = 412 MeV with parallel kinematics.) 





This work 


Jin et al. \U\ 


HCA 


DWEEPYl ^ 


DWEEPY2^ 


First Maximum 


1.08 


1.08 


0.98 


1.92 


1.21 


Second Maximum 


1.14 


1.13 


1.03 







^Extracted from Ref. ||10| 



TABLE II. Comparison of our spectroscopic factors obtained with the ccl and cc2 operators 
and with other relativistic and nonrelativistic analyses. 







ccl 


cc2 


Jin et al. 


121 


McDermott 


111 


nonrel. 


208 p]^ 


3Sl/2 
2(^3/2 


0.65(4) 
0.66(4) 


0.70(4) 
0.73(4) 


0.71 




0.65 




0.50(5) 
0.53(4) 


40Ca 


2Sl/2 
l'^3/2 


0.44(3) 
0.60(3) 


0.51(3) 
0.76(4) 


0.75 
0.80 




0.66 




0.50(6) ^ 
0.65(5) ^ 



^From Ref. g 
^Fvom Ref. IfTtl 
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